The Larmor clock and anomalous spin dephasing in silicon 
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Drift-difFusion theory - which fully describes charge transport in semiconductors - is also uni- 
versally used to model transport of spin-polarized electrons in the presence of longitudinal electric 
fields. By transforming spin transit time into spin orientation with precession (a technique called 
the "Larmor clock") in current-sensing vertical-transport intrinsic Si devices, we show that spin dif- 
fusion (and concomitant spin dephasing) can be greatly enhanced with respect to charge diffusion, 
in direct contrast to predictions of spin Coulomb-drag diffusion suppression. 

PACS numbers: 72.25.Dc 



Observation of spin precession and dephasing is not 
only clear evidence for noneouilibrium spin injection 
into nonmagnetic materials [l|, y], but it also provides 
a means to quantitatively measure electron transport 
dynamics. [3| As shown recently, spin precession measure- 
ments in a perpendicular magnetic field can be used to 
directly derive the empirical spin transit time distribu- 
tion via the Fourier transform, y This ability to extract 
time-of-flight data from a quasistatic measurement is en- 
abled by the encoding of spin transit time t into spin 
precession angle 9 = cot (where lo is spin precession fre- 
quency), a technique called the "Larmor clock" in the 
context of transport theory [5|-l8|. 

The power of this method to reveal the details of spin 
transport is diminished, however, by systematic sources 
of spin "dephasing" (uncertainty in spin orientation, A£). 
In lateral voltage-sensing|9l-ll2j or current-sensing[J, |l3| 
devices, for example, there is strong geometric spin de- 
phasing due to transit-length uncertainty induced by 
nonzero injector and detector width along the transport 
channel dimension. Vertical-transport devices, on the 
other hand, have true transverse symmetry, eliminating 
this systematic geometric spin dephasing and allowing 
the transport distribution to expose the consequences of 
the fundamental transport phenomena unencumbered by 
spurious effects lJ-[l6J . 

In this Letter, we use the concept of the Larmor clock 
to identify a phenomenon of anomalous spin diffusion, 
which behaves in a manner strikingly different from both 
the expectations of charge transport in the nondegener- 
ate regime and our present understanding of spin trans- 
port in semiconductors. Several potential physical origins 
are discussed, highlighting the need for a more sophisti- 
cated drift-diffusion theory necessary to accurately model 
nonequilibrium spin transport in nonmagnetic materials. 

Fig. [1] shows an example of spin current transit-time 
distribution extraction at a fixed temperature T =85K 
for different internal electric fields E=180 - 580 V/cm in 
vertical-transport L =350/im-thick Si(lOO) devices utiliz- 
ing ballistic hot electron injection and dctection[lJ, |l5|, 
Il7| (For a review of the device details and related tech- 
niques, see Ref. Il8|.) The inset shows raw spin pre- 
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FIG. 1: Transit-time distributions for spin current traveling 
through 350/im undoped Si under different internal electric 
fields (180-580 V/cm) at a temperature of 85K. Inset: Sym- 
metrized spin precession data, from which these distributions 
were calculated using the Fourier transform method described 
in Ref. |1 



cession data, symmetrized to avoid the effects of ferro- 
magnetic contact switching and maintain a real-valued 
Fourier transform. We control the spin precession fre- 
quency through the transverse magnetic field B, via 
uj — gfj.BB/h, where g is the electron g-factor, fiB is the 
Bohr magneton, and h is the reduced Planck constant. 
Oscillations as a function of magnetic field with average 
period 2Trh/gfiBt, where i is average transit time 17|, in- 
dicate coherent spin precession. The amplitude suppres- 
sion of these oscillations at higher magnetic field is the 
direct result of spin dephasing from nonzero A9 due to 
spin diffusion. As expected from basic charge transport, 
higher drift electric fields result in shorter average transit 
times and lower transit time uncertainty. These measure- 
ments can be repeated at different temperatures, and the 
normalized distributions obtained {P{t)) can be analyzed 
to determine any distribution moment. 

Most important here is the transit-time un- 



certainty (distribution standard deviation At — 
J j^ P{t)(t — tYdt), which is a direct measure of spin 
dephasing through A0 = wAi. This quantity is plotted 
in Fig. [21 along with similar data measured at other 
temperatures between 20K and 150K. All constant- 
temperature data show a similar trend of increasing 
uncertainty as the mean increases in smaller drift fields. 
The model universally used to simulate spin precession 
measurements like these incorporates only drift (due to 
an electric field E) at velocity u, diffusion (caused by ran- 
dom thermal motion) with coefficient D, and spin relax- 
ation with timescale r. The resulting Fokker-Planck rate 
equation is nominally identical to the minority-carrier 
equation used to successfully model charge transport in 
semiconductors, [13] except that it models the density of 
a single component of the spin vector rather than charge 
density. Due to the absorbing boundary conditions rele- 
vant for our current-sensing devices, the transit-time dis- 
tribution determined from the Green's function s(x, t) of 
this scalar drift-diffusion- relaxation equation is[20| 
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FIG. 2: Comparison of experimental and theory spin current 
transit time uncertainty (standard deviation of transit time 
distributions). Larger uncertainty corresponds to higher spin 
dephasing. Theory results rely on the drift-diffusion Green's 
function (Eqn. [l]), with Canali et al's data in Ref. [22 and the 
nondegenerate Einstein relation for accurate diffusion coeffi- 
cient modeling. 
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We expect our empirically measured distributions to 
have this form, for some value of D, v, and r. Spin life- 
time T has been measured in undoped Si 15| and its tem- 
perature dependence has been analyzed theoretically 2l| ; 
typical values are much longer than the transport times 
measured (< 60ns in Fig. [2]) and so have negligible effect 
on the transit-time uncertainty. Because electron density 
is in the nondegenerate regime in our experiments, it may 
appear reasonable to apply the fluctuation-dissipation 
theorem (Einstein relation, D — ^kBT/q, where ks is 
the Boltzmann constant and q is the fundamental elec- 
tron charge) which provides a relationship between D 
and V via the mobility /i — v/E. 

Using this simple theory, we can generate predictions 
of spin dephasing at constant D, as shown by the dot- 
dash and dashed lines (corresponding to 50cm^/s and 
500cm^/s, respectively) in Fig. [2 It is clear that the em- 
pirical data shown (labeled "experiment" ) have effective 
diffusion coefficients in excess of the latter, larger, value. 
When the Einstein relation (along with well-established 
charge transport data from Ref. [22l |) is used, uncertain- 
ties corresponding to diffusion coefficients of much lower 
magnitude (labeled "theory" in Fig. [5]) are predicted. 

What is the cause of this experimentally-observed en- 
hancement of spin dephasing? It has been known for 
some time that spin Coulomb drag[23|, |2J] causes an in- 
equality between charge diffusion and spin diffusion co- 
efficients and inapplicability of the Einstein relation, but 
always suppressing spin diffusion with respect to charge 
diffusion. Here, however, we see an apparently large 
enhancement. Because the Larmor clock depends on 



a well-defined precession frequency w, uncertainties in 
this quantity Aw, caused by inhomogeneous magnetic 
fields or g-factor variations, can certainly contribute to 
dephasing. However, we have confirmed that the drastic 
spin diffusion enhancement seen here would require inho- 
mogeneities of approximately 50%, which is inconsistent 
with our experimental conditions. 

Even more striking evidence for anomalous dephasing 
is presented in Fig. [3l where we show temperature de- 
pendence at a fixed electric field of E=580 V/cm. The 
transit-time distributions are shown in Fig. [3Ka), along 
with the corresponding symmetrized spin precession data 
in the inset. The spin dephasing temperature dependence 
(distribution standard deviation) is shown in Fig. [31(b). 
Dephasing is small at the lowest temperature (20K), but 
rises to a maximum near 60K, after which it slowly falls 
as temperature increases toward lOOK. Even without 
calculating the standard deviation from corresponding 
transit-time distributions, it is clear from the raw data 
in the inset of Fig. ^a.) that the width of the envelope 
function of the precession oscillations - inversely propor- 
tional to the dephasing - has a minimum near 60K. This 
behavior cannot be explained within the drift-diffusion 
model using the fiuctuation-dissipation theorem, which 
predicts a monotonic increase in dephasing as a function 
of temperature. [16| This prediction is shown using both 
the empirical mobility (u = L/iE) and charge-transport 
data (again, from Ref. |22|) as dashed and dotted lines, 
respectively; they are nearly identical. Once again, in 
addition to the monotonic nature of this prediction, the 
magnitude is far below what is seen experimentally. Sys- 
tematic sources of dephasing are temperature indepen- 
dent^ so the phenomenon seen here is clearly intrinsic to 



the details of spin precession during transport in Si. 

This nonmonotonic behavior is somewhat suggestive of 
thermally-activated electron trapping during transport. 
In such a model, average electron dwell times (rt) before 
emission from shallow traps are exponentially dependent 
on temperature according to the Arrhenius law: at high 
temperature, the trap time is very short, so electron trap- 
ping has negligible effect on the transit time distribution 
compared with the diffusion-induced distribution width. 
At low temperatures, the trapping time can be much 
longer than the spin lifetime, so when a trapped elec- 
tron is eventually released, it has become depolarized and 
makes no contribution to the overall spin signal. How- 
ever, at intermediate temperatures, the effects of electron 
trapping on distribution broadening may be significant. 

Let us now evaluate the applicability of this model. 
Due to the unimodal form of the distributions at inter- 
mediate temperatures shown in Fig. Elja), nearly all of 
the electrons would have to be affected by such traps, 
so the expected number of average trapping events k is 
large. In this limit, the convolution of k exponential- 
decay trap escape distributions (the Erlang distribution 
G{t;k,Tt)) approaches a gaussian form with mean to- 
tal trapping time krt and standard deviation vkTt |25| . 
Therefore, if this trapping model is applicable, the de- 
viation of the mean spin transit time from the expected 
charge transport time should be more than the difference 
in the observed standard deviation from that predicted 
using charge transport data for any fc > 1. However, as 
shown by the mean transit time data (grey, right axis) in 
Fig. EJb), this is not the case; spin and charge mobility 
are nearly identical. 

Furthermore, as shown in Fig. |4l increasing the inter- 
nal electric field only a nominal amount from 580V/cm 
to 1150V/cm strongly suppresses this temperature- 
dependent anomalous spin dephasing far beyond what is 
expected from trap ionization. Shallow trapping is thus 
not the cause of our observations of enhanced spin diffu- 
sion. The observed anomalous dephasing may therefore 
be regarded as the result of a more fundamental, intrin- 
sic, cause. For instance, due to the direct relationship 
between electron density and transit time in our insulat- 
ing samples [16|, higher electric fields result in more dilute 
electron densities and modify the electron-electron inter- 
action which may play a role in spin dephasing. 

Because spin is a vector quantity, modeling transport 
as a scalar as is typically done in the drift-diffusion equa- 
tion is necessarily incomplete. A more complete, three- 
component vector drift-diffusion theory would account 
for this, including tensor representation of the spin dif- 
fusion coefficient. [28[ It remains to be seen, however, 
whether the anomalous behavior described here would be 
captured solely by this more sophisticated framework. 

In addition, the interaction of phonons with spin- 
polarized electrons may become relevant at the transport 
lengthscale investigated here. In light of recent theoret- 



a) 



ffl 
O 




■ 580V/cni -200 200 

I . I Perpendicular Magnetic Field [Oe1 



5 10 15 
Transit Time [ns] 




0.0 



580V/cnn 



drift-diftusion theory 



20 40 60 80 

Temperature [K] 



100 



FIG. 3: Temperature dependence of spin dephasing in un- 
doped Si showing (a) spin current transit time distributions 
obtained from spin precession measurements (shown in inset, 
where the curved arrow emphasizes the extent of the envelope 
function modulating the coherent oscillations) in an internal 
electric field of 580 V/cm. (b) Temperature dependence of 
transit time uncertainty (standard deviation of the distribu- 
tions) showing nonmonotonic behavior that contrasts with 
the monotonic increase expected from drift-diffusion theory 
(Eqn. [l]) using the Einstein relation. Long dashes use the 
empirical mobility, whereas short dashes use charge transport 
data from Ref. [22| ]. Grey data (and right axis) shows the 
mean transit time of the empirical distribution (open circles, 
"spin") and prediction from Ref. ^^ (solid line, "charge"). 



ical developments with sophisticated calculations of the 
effects of spin polarized electron-phonon interactions uti- 

it is suggestive that 



lizing detailed bandstructure|21 



thermal conductivity due to three-phonon Umklapp pro- 
cesses also show nonmonotonic behavior [30| with a peak 
near 60K as in Fig. EJb) . Although study of the Elliott- 
Yafet electron-phonon interaction has focused solely on 
spin relaxation [2 llls Ij . it may be that similar physics has 
an effect on spin dephasing when all three components of 
spin are taken into account in a transport scenario. 

Finally, we wish to point out that this anomalous de- 
phasing behavior may be a universal characteristic of 
nonequilibrium spin transport processes. A notable ex- 
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FIG. 4: Suppression of anomalous spin dephasing with elec- 
tric field from 580 V/cm to 1150V/cm. Inset: Drift-diffusion 
predictions using the Einstein relation. Solid lines use the 
empirical mobility, whereas dashed lines use charge transport 
data from Ref. i22il. 



ample is evident in optical measurements of the per- 
sistent spin helix in two-dimensional GaAs32], where 
nonmonotonic temperature behavior of the spin-orbit- 
enhanced mode lifetime is observed. Because the "per- 
sistence" of the spin helix relies on control of spin dif- 
fusion, the similarity in behavior with what we present 
here in Figs. [3] and |3] is suggestive of analogous under- 
lying physics despite very different materials and even 
dimensionality. 26|, |27 1 

In conclusion, we have shown how the Larmor clock 
concept can be used as a powerful probe of spin dynamics 
in Si by providing a means to encode fast time-dependent 
transport processes into spin orientation that is measured 
quasistatically. A drastic enhancement of spin dephasing 
seen especially in nonmonotonic temperature dependence 
(from devices where geometric dephasing is eliminated by 
transverse symmetry) cannot be explained by our present 
understanding of spin transport which would contrar- 
ily predict a suppression of spin dephasing due to spin 
Coulomb drag. The challenges to construction of a suc- 
cessful theory capturing the observed phenomenon have 
been discussed; it is hoped that the suggestions here pro- 
vide a theoretical starting point for a more sophisticated 
model of nonequilibrium spin transport in nonmagnetic 
semiconductors . 
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